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Summary 

Recently reported experimental results demonstrate that 
restabilization of the low-Reynolds-number flow past a 
circular cylinder can be achieved by the placement of a 
smaller cylinder in the wake of the first at particular 
locations. Traditional numerical procedures for modeling 
such phenomena are computationally expensive. An 
approach is presented here in which the properties of the 
adjoint solutions to the linearized equations of motion 
are exploited to map quickly the best positions for the 
small cylinder's placement. Comparisons with experi- 
ment and previous computations are favorable. The 
approach is shown to be applicable to general flows, 
illustrating how strongly control mechanisms that 
involve sources of momentum couple to unstable (or 
stable) modes of the system. 


Introduction 

We present here a linear analysis that permits, through 
use of stability theory, an approximation to the shift in 
the stability characteristics of steady, two-dimensional 
(2D) flow past a circular cylinder that occurs as a result 
of the presence of an additional small cylinder within the 
flow. The flow under study is incompressible and 
viscous, but the analysis could undoubtedly be 
generalized to include compressibility as well as other 
effects. For this analysis the feedback mechanism is 
provided by the small cylinder, which in effect adds 
momentum to the fluid. Other types of control laws 
based upon unsteady momentum addition could be 
modeled using the kind of analysis presented here and 
applied to three-dimensional (3D) flows. 

The global dynamics of a system can be modified by 
specifying external sources which act in a way that 
depends upon the state of the system. The natural modes 
of the new system will differ from those of the old. If 
the system is not too highly distorted by the modifica- 
tion, we might expect to find that the natural frequencies 
of the new system have not been changed greatly. For 
linear waves in a continuous medium, adjoint solutions 
provide a means of examining in a simple fashion the 
effect of a particular feedback mechanism upon a particu- 
lar mode of the system. Should the mode be unstable, 
then the adjoint solution has particular use in pinpointing 
a mechanism for its restabilization. The formulation of 
the adjoint linear problem is a classical method in differ- 
ential equations, but because it is central to the calcula- 
tion we take time here to outline its construction with 
special reference to the incompressible Navier-Stokes 
equations. 


Viscous, incompressible, 2D flow past a circular cylinder 
is unstable above a critical Reynolds number of about 45 
based upon cylinder diameter, and stable below the criti- 
cal condition. The problem solved here is for a Reynolds 
number of 50, and the regular and adjoint eigensolutions 
are found for the unstable mode using a modification of 
the spectral approach of Zebib (ref. 1 ). The unstable 
waves in this case are believed to grow up to form the 
time-periodic Karman vortex street. Strykowski and 
Sreenivasan (ref. 2) performed both experimental and 
numerical examinations of the effect of introducing a 
smaller cylinder at various points within the wake of a 
larger one past which the flow is unstable. Suitable posi- 
tioning of the smaller cylinder was shown by both exper- 
imental observation and direct numerical simulation to 
suppress the onset of the vortex street. We present a sim- 
pler means of modeling the effect of the small cylinder 
in the wake of the first and demonstrate the restabilizing 
and, indeed, destabilizing effect which the small cylinder 
can have, by use of the adjoint to the unstable mode. The 
results of reference 2 are essentially verified, and 
furthermore there is a suggestion that there may be addi- 
tional regions, overlooked in reference 2, where the small 
cylinder could be positioned to restabilize the flow. 
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Adjoint Eigensolutions 

Adjoint solutions are used in the study of partial differ- 
ential equations (ref. 3), and find particular use in bifur- 
cation theory (ref. 4). In this section the relationships 
between linearized solutions and their adjoints are 
outlined within the framework of the incompressible 
Navier-Stokes equations, and their basic mathematical 
properties are presented. The results in this section are 
not mathematically new, but are included for 
completeness. 


Let us consider the steady, 2D, incompressible viscous 
flow field past a circular cylinder. (The application of 
this approach to time-periodic flows is covered in refer- 
ence 5.) The velocity field, V(r), is normalized by the 
uniform flow at infinity, Voo, and pressure, P(r), normal- 
ized by pVoo 2 , where p is the fluid density. If V\ is the 
component of V in the ith coordinate direction, then we 
have for the equations of motion 
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where R is the Reynolds number based upon cylinder 
diameter. Lengths are scaled on the cylinder radius. 
Numerically the problem is solved for the stream func- 
tion on an annular domain centered on the cylinder center, 
using a spectral method as described in reference 1. The 
reader is referred to this paper for details of the 
computational method. 

The flow shown in figure 1 is the stream function, ¥, for 
a Reynolds number of 50, and has been constructed using 
44 Chebychev polynomials in the radial direction and 
40 Fourier components in the angular direction. The outer 
computational perimeter is at 20 cylinder radii. The 
steady velocity field is then V = V x (Tz), where z is a 
unit vector normal to the plane of the 2D flow. The sepa- 
ration angle and the wake length agree with the results of 
reference 6. 

If we superimpose linear time-dependent disturbances 
upon this field, then they satisfy 
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where vj is the ith component of the velocity disturb- 
ance, v, and p is the pressure disturbance. Here the time 
variable is scaled by r/Voo, where r is the cylinder radius. 
The ith component of the linear operator L(V ; R) is 
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Taking the Fourier transforms in time, 
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gives 
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and 
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Applying boundary conditions v\ = 0 on the boundary of 
the computational domain defines an eigenvalue problem. 
(In our problem the conditions are applied on the 
cylinder surface and at 20 cylinder radii.) 

We will now formulate the adjoint equations. Central to 
the creation of an adjoint solution is the definition of an 
inner product. For a pair of complex vector fields, u and 
v, defined over the flow domain, we define the inner 
product to be 

[u,v] = j*u vdV (9) 

v 

after reference 7. We now look for an operator L*(V; R) 
for which 

[L(V;R)u,v] = [u,L*(Y;R)y] (10) 


for all u and v, where u satisfies the boundary 
conditions for the regular linear problem and v has 
boundary conditions chosen to ensure the validity of 
equation (10), 

From equations (7) and (8), the linear adjoint problem 
with respect to the inner product (eq. (9)) can now be 
formulated as 
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and 


3^ = 0 (12) 

With the boundary conditions chosen to ensure equa- 
tion (10), the adjoint solutions can be solved to give a set 
of eigenfrequencies and eigensolutions. If y( n ) is the nth 
solution to equations (7) and (8) with eigenfrequency 
co n , and y( m )* is the mth adjoint eigensolution to equa- 
tions (11) and (12) with frequency co^ , then it is easily 
shown that 

|co n ,y (m) *] = 0 (13) 

Consequently either the frequencies of the nth and mth 
mode are the same, or the regular and adjoint solutions 
are orthogonal with respect to the inner product. We can 
relabel the frequencies, order them in increasing order of 
stability (-lm(coo) > -Im(O)i) >...), and normalize the 
velocity fields so that 

co n = .and [O n > , 0 m >*] = 8 n m 

for n.m = 0,1,2,. . . (14) 


(8) 



The adjoint linear operator L*(V ; R) to the linearized 
Navier-Stokes equations operator L(V ; R) can be 
developed from 
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The Control of Unstable Flows 

The recent paper by Strykowski and Sreenivasan (ref. 2) 
documents the effect on overall system stability of plac- 
ing a small cylinder in the wake of a larger cylinder past 
which the flow is unstable. We model the effect of the 
small cylinder as a pointwise supply of momentum to 
the flow equal and opposite to the anticipated drag. A 
quasi-steady approximation to the drag upon a small 
cylinder in a slowly varying flow is used, and is believed 
to be valid since the time scale of adjustment of the flow 
behind the small cylinder is shorter than that behind the 
large cylinder by at least a factor of the diameter ratio. 

If V + v is the flow velocity in the vicinity of the small 
cylinder, then the drag vector is taken to be 


and 


d = e(V + v) | V + v| (19) 
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Then 


[L(V;R)u,y] = [u,L*(V;R)v] + Jfjdnj (18) 


A no-slip condition on the cylinder surface and at the 
outer computational boundary for both the regular and 
adjoint velocity fields ensures that fj = 0. 

At a Reynolds number of 50, the incompressible flow 
past a circular cylinder is just unstable, with the distur- 
bance being asymmetric and the eigenfrequency being 
coq = 0.376 - i 0.000353. The Strouhal number, based on 
cylinder diameter and with frequency f in Hertz, is then 
2fr/Voo = Re(o>o/7t) = 0. 120, as expected. The real and 
imaginary parts of the stream function are shown in fig- 
ures 2(a) and 2(b). Here 34 Chebychev polynomials have 
been used in the radial direction, and 35 Fourier compo- 
nents in the angular direction. (These have been calculated 
using a modification of the method of reference 1 which 
is not described here.) This clearly resembles the forms 
found in reference 8. Figures 3(a) and 3(b) show the 
corresponding adjoint solutions using the same number 
of coefficients. Most notable here is the occurrence of 
features upstream of the cylinder (the steady flow is 
from left to right), something not seen in the regular 
solutions. The adjoint solution is largest in the 
immediate wake of the cylinder. 


where £ represents pCod/2 for the small cylinder, 
where d is its diameter and C& is the effective drag coef- 
ficient. We assume £ is small. Taking the Fourier trans- 
form in time of the drag, and expanding in powers of the 
velocity disturbance, we find that 
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We model the presence of the small cylinder in the wake 
of the larger cylinder at position £q as a momentum 
supply, which has both a steady and unsteady component, 
on the right-hand side of the momentum equations of 
strength 


-e 


I V|V8(0)) + t=tV- v + | V|v 8(r-r 0 ) (21 


If we write V^) tor the steady flow past the cylinder, 
satisfying equations (1) and (2), then the presence of the 
small cylinder will lead to a modification of the steady 
flow component. If we write the new steady flow as 

V= V< 0) + eV (1) +0(e 2 ) (22) 


then it is not hard to show that we must have 

l(v(°>;R)v ( 1 ) =-| V (0) |y (0 >6(r-r 0 ) (23) 

The linearized equations of unsteady motion around this 
new base flow are 
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The value of EG)’ is the first-order correction to the 
eigenfrequency COq resulting from the placement of a 
small cylinder at position r = ro. The shift in the damp- 
ing rate is given by -£lm(co'); a negative shift encourages 
restabilization, whereas a positive shift enhances the 
instability. 


The choice of equation (21) for the momentum supply is 
specific to our present problem. Clearly this could repre- 
sent other kinds of control mechanisms in which momen- 
tum is supplied to the system. 

It is anticipated that there will be an eigenvalue of the 
inhomogeneous momentum equation (eq. (24)) that is 
close to the unstable eigenvalue for the flow past the 
large cylinder on its own. Consequently we write 

oo 

v = + b n v^ 4- O(e^) (25) 

n=l 

and 

CO = COq + £C0’ +0(£ 2 ) (26) 

where is the unstable eigensolution with eigen- 
frequency coo, for some constants b n and co\ Substituting 
into the inhomogeneous momentum equation and collect- 
ing terms of equal order in £, to zeroth order, we retrieve 
the unmodified eigensolution as expected. To order £, 

oo oo 
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From this point onward we ignore the influence upon the 
eigenvalue of the correction, to the steady flow and 
set A equal to zero. Although this factor may not be 
insignificant, it is ignored here for the sake of simplicity. 
The result appears to justify this simplification, but the 
point requires further investigation. 

The terms -Im(co') and Re(co’) are plotted in figures 4 
and 5, respectively, (with A = 0). The demarcation of the 
stability boundary can be achieved only with some 
determination of £. We do not attempt this, and restrict 
our interest purely to the relative merits of one position 
for the small cylinder versus another. Regions of large 
negative shift imply that placing the small cylinder in 
this area will have a strong stabilizing effect upon the 
flow. It is interesting to make the comparison with the 
hard-earned numerical results of reference 2; note that 
essentially the same features have appeared in our plot. 
Furthermore, an additional region where vortex shedding 
may be suppressed is indicated 7 to 9 cylinder diameters 
downstream. This was not reported in reference 2. It is 
not entirely clear whether this is a physical effect or an 
artifact of the finite computational domain. The experi- 
mental results of reference 2 are for Reynolds numbers 
of order 80 and above, but show essentially the same 
features as the more negative contours of figure 4. 
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Multiplying the equation by the adjoint solution to v^ 0 ^ 
and integrating over space eliminates the summation 
factors because of the properties of the adjoints, revealing 
that 
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y(0) 

P) 


\ 

\K°) • +| |v<°> 

/ 


♦y< 0) * +A(28) 


where the term on the right-hand side in brackets is 
evaluated at the small-cylinder location r 0 ,and 


A further feature that can be seen in figure 4 is the region 
adjacent to the main cylinder walls in which the insta- 
bility is enhanced by the small cylinder. This is to be 
expected, because when the small cylinder is close to the 
large cylinder the effective Reynolds number is increased, 
thus encouraging the instability. 

The map of the shift in the Strouhal number caused by the 
small cylinder shows regions where the frequency is 
increased, whereas in other regions it is decreased. It was 
observed in reference 2 that the vortex shedding eventu- 
ally returns at higher Reynolds numbers, with the 
Strouhal number significantly reduced. In the regions of 
strong suppression indicated in figure 4 (contour levels 
of about -0.02), the corresponding frequency shift in fig- 
ure 5 is negative. Therefore the results here are consistent 
with this frequency-reduction trend. 
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Conclusions 

It has been shown that an approximate linear stability 
analysis can be employed to illustrate how the stability 
of low-Reynolds-number flow past a circular cylinder is 
modified by the placement of a smaller cylinder in its 
wake. The results are shown to compare favorably with 
the more direct numerical computation of reference 2. 
While the choice of “control mechanism” was mechani- 
cal in this configuration, linear stability analysis could 
be applied to more general types of control in which a 
fluid system receives an external supply of momentum 
according to some prescribed control law. 
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Figure 1. The steady stream function for viscous, incompressible flow past a circular cylinder at a Reynolds number 
of 50. 
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(b) 

Figure 2. The (a) real and (b) imaginary parts of the stream function representing the unstable eigensolution for viscous, 
incompressible flow past a circular cylinder at a Reynolds number of 50. 
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(b) 

Figure 3 . The (a) real and (b) imaginary parts of the adjoint stream function representing the unstable mode for flow 
past a circular cylinder at a Reynolds number of 50, 
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Figure 4. The shift in the damping rate , -lm(co'), as a function of the placement position for the small cylinder at a 
Reynolds number of 50. The chain dash line shows the computed results of reference 2 for cylinder placement for 
neutral stability at a Reynolds number of 48 (diameter of small cylinder is 1/10 that of large). 



Figure 5. The shift in Re ((o') (proportional to the Strouhal number) as a function of the placement position for the small 
cylinder at a Reynolds number of 50. 
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